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Analysis of a multi-dimensional structure 


The invention pertains to a method for analysis of a multi-dimensional 
structure which includes a tubular structure from two-dimensional datasets. 

5 Such a method for analysis is known from the paper '3D reconstruction of 

coronary arterial tree to optimise angiographic visualisation \ by S. James Chen and John D. 
Canroll in IEEE Tr.Med.Imag.l9(2000)3 18-336 

10 The known method concerns reconstruction of a three-dimensional coronary 

arterial tree from two-dimensional datasets in the form of a sequence of x-ray angiograms. 
Such angiograms are made by injecting a patient to be examined with a contrast agent and 
passing x-rays from an x-ray source through the patient. The known method employs 
acquisition of two standard x-ray angiogram cine image sequences at respective orientations 

15 or views of the x-ray source relative to the patient. A pair of these cine images relating to 
corresponding phases of the heart is identified. This pair of identified corresponding cine 
images functions as the initial two-dimensional datasets. From each of the cine images of the 
corresponding pair, there are constructed vessel hierarchies that are represented by 
mathematical hierarchical digraph. Subsequently, a transformation characterising a spatial 

20 relationship between the views of the pair of corresponding cine images is constructed. This 
transformation is calculated as an optimum estimate of 3D bifurcation points given the set of 
2D bifurcation points in the hierarchical digraphs. The optimum estimate is obtained as the 
least squares errors of the 3D bifurcation points given sets of 2D bifurcation points. Notably, 
the known method is employed by an x-ray examination apparatus having a moveable gantry 

25 which carries the x-ray source and an image intensifier. The gantry angulations determine the 
orientation of the cine images. Inaccurate values for gantry angulations are recorded in the x- 
ray examination apparatus. The recorded values of the gantry angulations are employed as an 
initial estimate of the transformation that is calculated from the hierarchical bi-graphs. 
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An object of the invention is to provide a method of analysis of an multi- 
dimensional structure including a tubular structure which produces more accurate results for 
the local sizes of the tubular structure. 


This object is achieved according to the invention by a method for analysis of 
a multi-dimensional structure which includes a tubular structure from two-dimensional 
datasets for respective pre-determined projection directions, the method comprising the steps 
of 

1 0 - identifying at least one pair of corresponding initial projected centre points of the tubular 
structure in two respective initial two-dimensional datasets, 

- identifying at least one further projected centre point corresponding to said initial 
projected centre points in at least one further two-dimensional dataset, 

- derive a three-dimensional spatial position of the centre point of the tubular structure 
1 5 corresponding to said proj ected centre points, 

- identify projected edges of the tubular structure in said initial two-dimensional datasets 
and in said further two-dimensional dataset near the respective projected centre points 
and 

- derive a local size of the tubular structure at the three-dimensional spatial position of the 
20 centre point of the tubular structure from said projected edges and the predetermined 

projection directions. 


The invention is based on a set of two-dimensional datasets for pre- 
25 determined projection directions. Notably, these two-dimensional datasets are formed as 
digitised x-ray images for respective predetermined projection directions. The projection 
directions are pre-determined as the x-ray projection images are generated by an x-ray 
examination apparatus having pre-calibrated orientations of the x-ray source and x-ray 
detector. Such a pre-calibrated x-ray examination apparatus is described in the international 
30 application EP01/12743 (PHNL000587). 

According to the present invention there are identified corresponding initial 
projected centre points of the tubular structure in two initial two-dimensional datasets. These 
projected centre points are corresponding in that they pertain to the same centre point in the 
tubular structure in the multi-dimensional dataset. The corresponding initial projected centre 
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points are for example identified on the basis of features in the two-dimensional dataset that 
are well recognised. For example, the tubular structure may relate to a part of a patient's 
system of bloodvessels, notably the coronary arterial vessel tree. In the two-dimensional 
datasets representing projections of the system of bloodvessels, the initial projected centre 
5 points can be accurately identified on the basis of anatomical information in the projection 
images forming the two-dimensional datasets. Notably, suitable anatomical features are 
bifurcation points of the branched vessel structure which are often well recognised in the 
two-dimensional datasets when displayed on a display screen. 

Subsequently, the position in three-dimensional geometric space of the centre 
1 0 point in the tubular structure is calculated from the pre-determined projection directions of 
the initial two-dimensional datasets and the locahon of the initial projected centre points in 
these initial two-dimensional datasets. To obtain accurate values for the local size of the 
tubular structure in at least one further projected two-dimensional dataset at least one further 
projected centre point is identified. This further projected centre point also corresponds to the 
1 5 centre point in the tubular structure that corresponds to the initial projected centre points. 
That is, the initial projected centre points and further centre points in one or several further 
two-dimensional datasets all pertain to the same centre point in the tubular structure. In the 
initial two-dimensional datasets and in the further two-dimensional datasets projected edges 
of the tubular structure are identified locally at the projected centre points in the respective 
20 two-dimensional datasets. There are various ways to identify the projected edges e.g on the 
basis of local differences between brightness values of grey values or local gradients of 
brightness or grey values ion the two-dimensional datasets. From these projected edges and 
the pre-determined projection directions of the respective two-dimensional datasets an 
accurate estimate is computed of notably the cross-sectional size of the tubular structure. In 
25 particular the projected edges in the respective two-dimensional datasets and their pre- 
determined projection direction yield an inscribed polygon of the local cross-sectional area of 
the tubular structure locally at the centre point at issue. Consequently, the more two- 
dimensional datasets for different projection directions are employed, the larger the number 
of vertices of the inscribed polygon and the more accurate the estimate of the local cross 
30 section of the tubular structure on the basis of the inscribed polygon. 


These and other aspects of the invention will be further elaborated with 
reference to the embodiments defined in the dependent Claims. 
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Preferably on the basis of the projection directions of the initial two- 
dimensional datasets and the spatial position of the centre point the first epipolar line is 
computed as the intersection of the projection plane of the second initial two-dimensional 
dataset and the plane spanned by the projection directions of the initial two-dimensional 

5 dataset and containing the centre point at issue in the tubular structure. The initial projected 
centre point in the second initial two-dimensional dataset is easily identified on the epipolar 
line, since only positions on the epipolar line need be considered. Moreover, the position of 
the projected centre point on the epipolar line is obtained very accurately. The accuracy of 
the position of the projected centre point in the second initial two-dimensional dataset is 

1 0 determined by the accuracy of the calibration of the pre-determined projection directions. It 
appears that the pre-determined projection directions can be calibrated at such accuracy that 
the position of the projected centre point is determined far more accurately than the pixel-size 
in the two-dimensional dataset. 

It is noted that further epipolar lines can be computed for any further pair of 

15 two-dimensional dataset for respective projection directions. Projected centre points are 

located on such epipolar lines so that identification is easily done since only positions on the 
epipolar line need be considered. Projected centre points can also be identified as 
intersections of epipolar lines pertaining to respective pairs of projection directions. 

In a further preferred implementation of the method according to the 

20 invention, projected centre lines of the projections of the tubular structure are constructed. 
These projected centre lines represent the projections of the local longitudinal axes of the 
tubular structure. For example the projected centre lines represent the central axis along the 
lumen of the blood vessel at issue as projected in the x-ray projection image. The projected 
centre lines are obtained from respective series of projected centre points in the individual 

25 two-dimensional datasets. Such series of projected centre points can be obtained in the same 
way as discussed above for the corresponding initial projected centre point and further 
projected centre points. The projected centre lines are for example found by connecting 
subsequent projected centre points in the individual two-dimensional datasets. hi individual 
two-dimensional datasets the projected centre lines are corrected on the basis of the 

30 datavalues in the two-dimensional dataset at issue. For example the correction of the 

projected centre lines may be carried-out on the basis of a distribution of datavalues near the 
projected centre lines, on the basis of local gradients of the datavalues, local differences 
between datavalues or on the local Hessian matrix (i.e. local spatial 2 nd derivatives) of the 
datavalues near the projected centre line. Thus, on the basis of the datavalues close to the 
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projected centre line, the location of the centre line is corrected. Consequently, the positions 
of the projected centre points are corrected and made more accurate. Thus, also the accuracy 
of the local size of the tubular structure is improved. 

In particular the local size of the tubular structure is calculated as the local 

5 diameter at the centre point at issue and transversely to the corrected projected centre line. In 
this way the local size of the tubular structure is calculated in the respective projection 
directions. In this way complicated shapes, different from a circular form, of the cross section 
of the tubular structure is taken into account 

Accurate results for in particular the projected centre line(s) are found when 

10 the initial two-dimensional datasets have projection directions that enclose an angle in the 
range of 30° to 150°. Very accurate results are found when the projection directions of the 
initial two-dimensional datasets differ about 90°. Preferably, the projection directions of the 
further two-dimensional datasets enclose an angle also within the range between 30° to 150°. 
Consequently, the range of projection directions for which the two-dimensional dataset need 

15 to be acquired is less than the full range of 1 80° so that accordingly the time to acquire the 
two-dimensional datasets is reduced. Further, it appears that the angle enclosed between the 
projection directions of the initial two-dimensional datasets may be selected anywhere in the 
range of 30° to 1 50°, an in particular cumbersome setting said angle at 90° is not necessary. 
However, it is noted for projections that differ about 90°, overlap and foreshortening of the 

20 projected tubular structure are relatively less so that corresponding centre points are easily 
identified. 

Preferably, in particular when (almost) periodic motion occurs in the tubular 
structure time gated acquisition of the two-dimensional dataset is employed. In time gated 
acquisition the datavalues of the respective two-dimensional datasets are acquired in 

25 successive time periods in which the state of motion of the tubular structure is essentially the 
same. For example successive two-dimensional datasets at respective projection directions at 
periodic instants in time at which essentially the same motion state recurs. Consequently, the 
two-dimensional datasets form in effect a stroboscopic representation of the tubular structure 
at respective projection directions. For example when the tubular structure pertains to the 

30 system of bloodvessels of the patient's heart, time gated acquisition can be carried out by 
ECG-triggered acquisition. 

The method of the invention, in particular when applied to the bloodvessel 
system of the patient to be examined leads to an adequate spatial rendition of the bloodvessel 
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structure which is a useful technical tool for the physician to examine the patient's 
bloodvessel system. 

The invention further relates to a computer programme as defined in Claim 7. 
The computer programme according to the invention can be loaded in the working memory 

5 of a computer so as to allow the computer to carry-out the method of the invention. The 
computer programme can be supplied on a data carrier such as a CD-ROM or the computer 
programme can be downloaded from a data network such as the world-wide web. The 
invention also relates to a workstation as defined in Claim 8. The workstation is arranged to 
carry-out the method of the invention and to that end the workstation is provided with 

10 hardware circuitry of a processor for cairying out the method steps or a computer programme 
according to the invention is stored in the working memory of the processor. Preferably, the 
workstation according to the invention is set up with an x-ray examination apparatus such 
that image data from the x-ray examination apparatus are supplied to the workstation as the 
two-dimensional datasets. 

1 5 These and other aspects of the invention will be elucidated with reference to 

the embodiments described hereinafter and with reference to the accompanying drawing 
wherein 


20 Figure 1 shows an x-ray examination apparatus in which the invention is 

employed, 

Figure 2 shows a diagrammatic representation of the identification of centre 
points in respective two-dimensional datasets on the basis of epipolar lines as employed in 
the method of the invention and 
25 Figure 3 shows a diagrammatic representation of the identification of 

projected edges of the tubular structure and of the derivation of the local size of the tubular 
structure as employed in the method of the invention. 

30 Figure 1 shows an x-ray examination apparatus in which the invention is 

employed. The x-ray examination apparatus comprises an x-ray source 2 and an x-ray 
detector 3 which are mounted on an inner C-arm 12. The inner C-artn 12 is mounted on a 
outer C-arm 14, such that the inner C-arm 12 can rotate around an axis 24 which is in the 
plane of the inner C-arm. The axis Z4 is often termed the *propellor-axis\ The outer C-arm 14 
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is mounted in a sleeve 15 which is suspended from the ceiling. The sleeve 15 is rotatable 
around axis Zi, which extends transverse to the long axis of the sleeve 15, in this example the 
axis zi is vertical. The outer C-ann 14 is moveable through the sleeve 15 so as to rotate the x- 
ray source 2 with the x-ray detector 3 around an axis Z3, consequently the propellor axis Z4 

5 can be set at an angle a to the horizontal axis z&. 

The object to be examined, such as a patient to be radiologically examined is 
placed on the patient table 4. Several two-dimensional datasets are made by positioning the 
x-ray source and x-ray detector according to respective predetermined projection directions. 
The pre-determined projection directions are achieved by displacing the outer C-arm 14 so as 

10 to set the propellor axis Z4 and to rotate the inner C-ann 12 around the propellor axis. 

The two-dimensional datasets are produced in the form of x-ray projection 
images, which are formed on the input face of the x-ray detector 3. Namely, owing to 
differences in the x-ray attenuation in the patient to be examined, spatial intensity variations 
occur in the x-ray beam having passed through the patient to be examined . These spatial 

1 5 brightness variations give rise to a different brightness values in the x-ray image formed on 
the input face of the x-ray detector 3. The x-ray detector may be an x-ray image intensifier 
having an input screen as the input face, or the x-ray detector may be a flat x-ray detector 
which includes a sensor matrix of x-ray sensitive sensor elements which form the input face. 
Preferably, the x-ray images are digitised in an analogue-to-digital converter that is provided 

20 in the output circuit of the x-ray detector 13. These digitised x-ray images form the two- 
dimensional datasets. 

Figure 2 shows a diagrammatic representation of the identification of centre 
points in respective two-dimensional datasets on the basis of epipolar lines as employed in 
the method of the invention. The two initial two-dimensional datasets 21,22 represent x-ray 

25 images formed by x-irradiating the patient's heart to be examined from respective positions/ 4 
and/ the x-ray source 2, in particular of the x-ray focus in the x-ray source. Accordingly, the 
x-ray images show respective projections of the patients coronary vessel structure. The 
contrast of the vessels is enhanced as an x-ray absorbing contrast agent has been administered 
to the patient's blood. Corresponding centre points C A and C B are indicted in the respective 

30 two-dimensional datasets 21,22. Notably, the corresponding centre points C A and Cb are 
chosen as an anatomically easily recognised as a bifurcation point in the vessel structure. 

At least one further two-dimensional dataset in the form of an x-ray image 3 
taken from x-ray source position^. The first and second epipolar lines E* and E? are shown. 
The first epipolar line E A is the line of intersection of the plane onto which the x-ray image 
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associated with two-dimensional dataset 23 is projected in the plane spanned by the 
projection directions of the first initial two-dimensional dataset 1 and of the further two- 
dimensional dataset 23. Similarly, the second epipolar line E A is the line of intersection of the 
plane onto which the x-ray image associated with two-dimensional dataset 23 is projected in 
5 the plane spanned by the projection directions of the second initial two-dimensional dataset 2 
and of the further two-dimensional dataset 23. As the centre points Ca and Cb both 
correspond to the same centre point in the patient's vessel structure, the corresponding centre 
point in the further two-dimensional dataset 23 is found as the intersection of both epipolar 
lines. 

10 Figure 3 shows a diagrammatic representation of the identification of 

projected edges of the tubular structure and of the derivation of the local size of the tubular 
structure as employed in the method of the invention. For individual elements of a vessel, the 
3D centreline points are calculated having a corresponding 2-D centreline point in every 
projection. For each of these points, two edge point pairs (Er, Ei) exist from which two radius 

1 5 vectors (Rr, Rj) can be reconstructed as follows. First the minimum distance between the line 
/ connecting the image point P B 2d and its corresponding focal spot./ 8 and the 3D centreline 
point P°3o is calculated. The corresponding point on the line L is defined as P^d- If the 3D 
centrelines calculation involves more than two projections the distance will not be zero due to 
the errors mentioned above. The two vectors pointing perpendicular from P B 3D to the edges of 

20 the triangle defined by the points Er, Ei and the focal spot f are the sought radius vectors Rr 
and Rj. The vessel lumen can be estimated from a single view by use of the measured vessel 
size in the image assuming circularity of the vessel lumen. When two views are acquired 
elliptical cross section can be employed. For more than two projections, interpolation is 
performed. The quality of this estimation increases with the number of projections used. 
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1 . A method for analysis of a multi-dimensional structure which includes a 
tubular structure from two-dimensional datasets for respective pre-determined projection 
directions, the method comprising the steps of 

- identifying at least one pair of corresponding initial projected centre points of the tubular 
5 structure in two respective initial two-dimensional datasets, 

- identifying at least one further projected centre point corresponding to said initial 
projected centre points in at least one further two-dimensional dataset, 

- derive a three-dimensional spatial position of the centre point of the tubular structure 
corresponding to said projected centre points, 

10 - identify projected edges of the tubular structure in said initial two-dimensional datasets 
and in said further two-dimensional dataset near the respective projected centre points 
and 

- derive a local size of the tubular structure at the three-dimensional spatial position of the 
centre point of the tubular structure from said projected edges and the predetermined 

1 5 projection directions. 

2. A method for analysis of a multi-dimensional structure as claimed in Claim 1, 
wherein identifying the at least one pair or corresponding initial projected centre points 
involves 

20 - identifying a first initial projected centre point in a first initial two-dimensional dataset 

- computing a first epipolar line as the intersection of a projection plane of a second initial 
two-dimensional dataset and the plane spanned by the projection directions of said initial 
two-dimensional datasets and 

- identifying a second initial projected centre point in the second initial two-dimensional 
25 dataset on the first epipolar line in the second initial two-dimensional dataset 

3. A method as claimed in Claim 1, wherein 

- in at least one of the two-dimensional datasets a series of projected centre points is 
identified 
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• in said at least one two-dimensional dataset a projected centre line is derived from said 
series of projected centre points 

- said projected centre line are corrected on the basis of datavalues in the at least one two- 
dimensional dataset. 

5 

4. A method as claimed in Claim 3, wherein the local size of the tubular structure 

is derived from an intersection of the identified projected edges in said at least one two- 
dimensional dataset and a direction transverse to the corrected projected centre line 

10 5. A method as claimed in Claim 1 , wherein the initial two-dimensional dataset 

have projection directions enclosing an angle in the range of 30° to 150°. 

6. A method as claimed in Claim 1, wherein the two-dimensional datasets are 
acquired by way of time-gated data-acquisition. 

15 

7. A computer programme comprising instructions for 

- identifying at least one pair of corresponding initial projected centre points of the tubular 
structure in two respective initial two-dimensional datasets, 

- identifying at least one further projected centre point corresponding to said initial 
20 projected centre points in at least one further two-dimensional dataset, 

- deriving a three-dimensional spatial position of the centre point of the tubular structure 
corresponding to said projected centre points, 

- identifying projected edges of the tubular structure in said initial two-dimensional 
datasets and in said further two-dimensional dataset near the respective projected centre 

25 points and 

- deriving a local size of the tubular structure at the three-dimensional spatial position of 
the centre point of the tubular structure from said projected edges and the predetermined 
projection directions. 

30 8. A workstation arranged to 

- identify at least one pair of corresponding initial projected centre points of the tubular 
structure in two respective initial two-dimensional datasets, 

- identify at least one further projected centre point corresponding to said initial projected 
centre points in at least one further two-dimensional dataset, 
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derive a three-dimensional spatial position of the centre point of the tubular structure 
corresponding to said projected centre points, 

Identify projected edges of the tubular structure in said initial two-dimensional datasets 
and in said further two-dimensional dataset near the respective projected centre points 
and 

derive a local size of the tubular structure at the three-dimensional spatial position of the 
centre point of the tubular structure from said projected edges and the predetermined 
projection directions. 
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